Quantum coherence and sensitivity of avian magnetoreception 
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Migratory birds and other species have the ability to navigate by sensing the geomagnetic field. 
Recent experiments indicate that the essential process in the navigation takes place in bird's eye 
and uses chemical reaction involving molecular ions with unpaired electron spins (radical pair). 
Sensing is achieved via geomagnetic-dependent dynamics of the spins of the unpaired electrons. 
Here we utilize the results of two behavioral experiments conducted on European Robins to argue 
that the average life-time of the radical pair is of the order of a microsecond and therefore agrees 
with experimental estimations of this parameter for cryptochrome — a pigment believed to form the 
radical pairs. We also found a reasonable parameter regime where sensitivity of the avian compass 
is enhanced by environmental noise, showing that long coherence time is not required for navigation 
and may even spoil it. 

PACS numbers: 03.67.-a,03.65.Yz,82.30.-b 



Recently there has been a growing interest in the ap- 
plication of quantum mechanics to understand many bi- 
ological phenomena such as photosynthesis \MB , process 
of olfaction 0, Q , enzymatic reactions fill Till] or avian 
magnetoreception [L2Tfl5j . These interests have brought 
physicists, chemists, and biologists at the same platform 
and led to the beginning of a new interdisciplinary sub- 
ject called quantum biology [HI [T3|. A major motiva- 
tion of these studies is to understand how nature utilizes 
purely quantum phenomena to optimize various biologi- 
cal processes. 

Here we are specifically interested in the avian mag- 
netoreception. It is very plausible that the navigation 
ability of some migratory birds is governed by the mecha- 
nism based on geomagnetic-dependent dynamics of spins 
of unpaired electrons in a radical pair. A recent theoreti- 
cal study has estimated both the life-time of the pair and 
the coherence time of this dynamics to be of the order of 
tens of microsecond [HI . The basic criterion used there 
postulates that bird's navigation is disturbed if the signal 
produced by the dynamics is independent of the orienta- 
tion of the geomagnetic field. This criterion together with 
the results of behavioral experiments in which European 
Robins could not navigate in a weak oscillating magnetic 
field [ill, [l9| led to the estimated life time and coher- 
ence time. Here we additionally take into account the re- 
sults of other behavioral experiments in which the same 
species were observed to be temporarily disoriented in 
a constant magnetic field sufficiently stronger or weaker 
than the geomagnetic field [2(1 Hl| . We estimate the life 
time and coherence time of the order of several microsec- 
onds. Our estimate is consistent with that obtained in 
a recent behavioral experiment [13l ] and also with the in 
vitro experiment using cryptochrome \22] , a pigment be- 
lieved to form the radical pairs. Furthermore, we demon- 
strate theoretically that environmental noise can enhance 
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the sensitivity of the avian compass, i.e. sensitivity in 
the presence of noise is better than without noise. This 
increase of sensitivity has resonant character and shows 
that long coherence times sometime may be disadvanta- 
geous for navigation. 

Let us begin with certain functional properties of avian 
compass which have been observed in different behavioral 
tests: (i) A very early experiment with European Robin 
showed that in contrast to the well-known magnetic ma- 
terial based physical compass the avian compass does not 
depend on the polarity of ma gne tic field but only on the 
inclination of magnetic field 23]; (ii) An experiment in 
which geomagnetic field was supplemented with a very 
weak radio-frequency (RF) field showed that birds dis- 
oriented at the frequency of the RF field which is reso- 
nant with the energetic splitting of a free electron (due 
to its spin) induced by the local geomagnetic field. This 
was observed only when the RF field was not parallel to 
the local geomagnetic field [l8|, [l9j] ; (iii) Another experi- 
ment showed that the avian compass works within a nar- 
row "functional window" around the local geomagnetic 
field. The compass ceased functioning when intensity of 
the magnetic field increased/decreased by about 30% of 
the local geomagnetic field. However, birds exposed to 
the new intensity of the magnetic field sufficiently long 
adapted themselves and their compass again worked cor- 
rectly mni. 

These experimental findings are consistent with the fol- 
lowing most widely used model of avian magnetorecep- 
tion (see Ref. 24 1 for a recent review) . In Ref. [25} , a new 
type of RP based model has very recently been proposed. 
It is assumed that bird's retina contains a photoreceptor 
pigment with molecular axis direction dependent on its 
position in the retina. Absorption of incident light by 
a part of the pigment results in electron transfer to a 
suitable nearby part and in this way a radical pair is 
formed, i.e. a pair of charged molecules each having an 
electron with unpaired spin. In the external magnetic 
field the state of electron spins undergoes singlet-triplet 
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transitions and at random times the radical pairs recom- 
bine forming singlet (triplet) chemical reaction products. 
The amount of these chemical products varies along the 
retina as the direction of molecular axis changes, and the 
shape of this profile is believed to be correlated in bird's 
brain with the orientation of the geomagnetic held. 

The mathematical model of this mechanism of a chem- 
ical reaction based avian compass involves coupling of 
electron spins to external magnetic held and to the spins 
of the molecular core. The property (ii) suggests that 
the spin of one of the electrons in the radical pair is ef- 
fectively uncoupled from spins of any other particles and 
only interacts with the external held. In a simple qual- 
itative model the spin of the other electron is coupled 
to effective spin-^ core of the molecular part ("nuclear" 
part) [15| . The corresponding Hamiltonian is as follows: 

H = T . A* . Sx+jt . (S 1 + S 2 ), (1) 

where / is the nuclear spin operator, A is the hyper- 
hne (HF) tensor, Si = {ax , Vy , &z} are the elec- 
tronic spin Pauli operators, it is the magnetic field, and 
7 = \p§g is the gyromagnetic ratio with /i being Bohr 
magneton and g the electronic g-factor. Here we assume 
that the <?-factors are the same for both electrons and 
set it at the value corresponding to the free electron, 
i.e. g = 2. Following Ref. flElj t he HF tensor is cho- 
sen axial and anisotropic, i.e., A = diagjA^, A y , A z } 
with A x = A y = a = A z /2. We define a quantity 
A = (A 2 + A^ + A 2 ) 1 / 2 as a measure of HF coupling 
strength and a = h/A as a measure of this strength 
in the time units, h stands for Planck constant. Here 
a = hj\f%a. 

In order to verify how this model recovers the proper- 
ties of behavioral experiments, one varies external mag- 
netic field 

it = B Q (sin 9 cos <fi, sin 9 sin <fi, cos 9) 

+ i? r f cos(wi) (sin a cos/3, sin a sin/3, cos a), 

where Bo gives the strength of the local geomagnetic field 
or the artificially changed constant field, and i? r f is the 
additional RF-field with frequency u>, which is switched 
on optionally. For all our numerics, the RF field is or- 
thogonal to the static field. Due to the axial symme- 
try of HF tensor we set <fi — without loss of gener- 
ality. We also assume that /3 = for the oscillating 
field's direction. Since most of the behavioral exper- 
iments were performed in Frankfurt, we set the mag- 
netic field strength \B \ = B = 47 pT as the refer- 
ence value. The corresponding Larmor precession period 
is h/2~/Bo = 0.76 /is.Q The HF interaction strength a 
typically acquires values in the range 10ns — 1 ps [24| . 



The factor of two emerging here is due to our definition of Si as 
Pauli operators. 



but we have verified that only for those radicals with 
the coupling strength close to the corresponding geomag- 
netic strength of 0.76 ps present model predicts working 
compass. Therefore, here we study two different cases 
of a = 1.0 ps and 0.5 ps, which are slightly greater and 
smaller than the geomagnetic strength. 

The dynamics of the system consisting of two electrons 
and one nuclei depends on two processes: coherent evolu- 
tion determined by the above Hamiltonian and recombi- 
nation of the pair to create spin-state dependent chemi- 
cal products. The spin-chemistry community models the 
above processes following the Haberkorn approach [26[ 
according to which the dynamics of the density matrix p 
describing the whole system is governed by equation 

^ = -i[H, p] - | (Q sP + P Q S ) - k - (Q T p + pQt), (3) 

where k determines the reaction rates for the singlet and 
triplet recombination, here assumed to be the same, and 
Q s and Q T are the projection operators onto the singlet 
and triplet subspaces. Although other treatments are 
possible [13, HH , here we follow the Haberkorn approach 
as it is the most consistent for the spin-selective recom- 
bination of radicals 29] . We begin the evolution at the 
moment of radical pair creation with the initial density 
matrix 

p(0) = ^lw®hTW"|. ( 4 ) 

where the electron pair is in the singlet state \ip~) — 
"^75 (I tl) ~ I It)): an d the nuclear spin is in a completely 

mixed state ^Iat. The singlet product yield is defined as 
the amount of product decaying via the singlet channel 

$s = k (il>-\Ti N p(t)\il>-)dt, (5) 
Jo 

where the partial trace of the density matrix is taken over 
the nucleus subspace. 

We proceed to estimate the parameter fc, inverse of 
which gives the average life-time of the radical pair. For 
this purpose, we use behavioral tests (ii) and (iii) which 
describe bird's disorientation under the influence of RF 
field and weaker artificial geomagnetic field. We shall not 
be interested here in the process of bird adaptation to the 
artificial magnetic field and only use the fact that the 
latter disturbs the compass. Nevertheless we would like 
to emphasize that adaptation is a physiological process, 
therefore its explanation does not warrant any modifi- 
cations of the present model. The model describes con- 
tribution to the chemical products from a single radical 
pair whereas the brain receives a cumulated signal from 
all created pairs (that are assumed to evolve indepen- 
dently). Under the change of magnetic field intensity the 
cells of the retina may start to vary the amount of cre- 
ated pairs and in this way control the intensity of the 
cumulated signal. Similar mechanism is present in verte- 
brate photoreceptors that can adapt their pigment sen- 
sitivity to external light by varying concentration of cal- 
cium ions [30] . As an experimental indication that these 
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FIG. 1. (Color online) Life-time estimation of a radical pair. 
Each plot presents angular dependency of the singlet yield 
in presence of local geomagnetic field of 47 pT (•), in 30% 
weaker field (■), and for the geomagnetic field augmented 
with the resonant RF field (a). Calculations are performed 
for two strengths of the hyperfine coupling: slightly smaller 
than the geomagnetic strength (a = 0.5 ps, three lower plots) 
and slight bigger than the geomagnetic strength (a = 1.0 ps, 
three upper plots). On plots from left to right we display the 
growing value of the recombination rate k used in numerics. 
The radical pair life-time is the inverse of k. Since birds are 
disoriented both in the weaker field and in the presence of 
the RF field we postulate that the criterion for estimating k 
should be based on the comparison of the profile of triangles 
to the profile of squares. This leads to the average life-time 
of the order of several microseconds (see main text for de- 
tails) in agreement with experiments on pigment believed to 
be present in the bird's retina. 



two adaptation phenomena may share common mecha- 
nism we note that the full avian compass adaptation to 
external magnetic field intensity takes about the same 
time as adaptation of vertebrate photoreceptors to ex- 
ternal light intensity. 

Fig. [1] presents the angular dependence of the singlet 
yield for different k values and for three different kinds 
of magnetic field environment: in the geomagnetic field 
of 47/zT, in 30% weaker magnetic field as in the experi- 
ments [2(], , and in the geomagnetic field augmented 
by weak resonant RF field of £? r f = 150nT as in the ex- 
periment [HI], and also for a fair comparison with Ref. 
15]. All the plots show that the singlet yield angular 
profile changes slightly when the magnetic field intensity 
is reduced by 30%. Note that there is a clear angular 
dependency of the profile both in the original and the 
weaker field. Since the behavioral tests showed that the 
weaker magnetic field disrupts the avian compass, the de- 
gree of flatness of the angular profile cannot be the only 
property responsible for avian navigation. The profile 
has to rather be within a proper range from the profile 
corresponding to the geomagnetic field alone. Therefore 
when the resonant RF is added to the geomagnetic field, 



it can disrupt the avian compass not only by completely 
washing out the angular dependency of the profile [15[ 
but also by changing the amount of singlet yield such 
that the profile is outside the range estimated by ap- 
plying the weaker magnetic field. We use this criterion 
to estimate the parameter k. Consider first three upper 
plots of Fig. [T]that present results of calculations for the 
HF strength of a = 1.0 ps. The left plot assumes k = 10 s 
and shows that the profile corresponding to the RF field 
is much below the profile corresponding to the weaker 
magnetic field. Similar behavior is observed for lower 
values of k which we have not presented here. In the 
middle plot (fc = 1.5 x 10 5 ) the RF field profile is slightly 
below the profile corresponding to the weaker field and 
therefore correctly predicts that bird's navigation is dis- 
turbed. The right plot of k = 2.5 x 10 5 shows to the 
contrary that the RF field profile is slightly above the 
profile corresponding to the weaker magnetic field, and 
therefore avian compass is not disturbed in contradiction 
with experiments. This implies that all values of k be- 
low 2 x 10 5 are consistent with experiments. We use the 
critical value k = 2 x 10 5 for further calculations. The 
corresponding average life-time of the radical pair is 5ps. 
Similar analysis applied to the three lower plots, for the 
case of a = 0.5 ps, gives estimation of k — 1.5 x 10 5 and 
the corresponding average life-time 6.7 ps. These esti- 
mations are come close to the experimental estimations 
for the cryptochrome [22l| . which is believed to be the 
photo-receptor pigment responsible for the radical pair 
based avian compass of European Robins; and they also 
come within the range 2ps — 10/xs, the lifetime estimated 
in a recent behavioral experiment [llj]. We have also 
done the same analysis using stronger (£? r f = 470nT) 
and weaker (B r f = A7nT) resonant RF-field, which are 
respectively 1% and 0.1% of the local geomagnetic field 
strength, and both have been reported to disturb avian 
compass [l3j. For these two cases, we estimate the life- 
time about 1 — 2ps and 11.0 — 12.5ps, respectively (3lj . 
These estimated lifetimes are within or very close to the 
above mentioned range. We also study the case where 
the geomagnetic field (47/iT) is replaced by the artifi- 
cially created stronger (94/iT) static field. The estimated 
life-time for this case is also within the above mentioned 
range [3l| . 

After setting the value of decay parameter k, we in- 
vestigate the effect of environment on the singlet yield. 
Following Ref. [l5| we describe the environment by a 
standard Lindblad formalism 32, 331: 



dp 
~db 



-i[H, p] - - (Q s p + pQ s ) - - (Q T p + pQrj 



LipL\ - - (l\l iP + P L\U 



(6) 



where there are six noise operators Li altogether, three 
acting on the first electron and three acting on the second 
electron. The three noise operators for each electron are 
taken to be the Pauli matrices a x , a y and a z . We assume 
that the noise parameters Ti are all equal and given by 
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FIG. 2. (Color online) Magnetic sensitivity of avian compass 
in presence of environmental noise. The upper plots present 
results of calculations for hyperfine coupling strength of a = 
1 [is and k — 2 x 10 5 , the lower plots are for a = 0.5 /zs and 
k — 1.5 x 10 s . The left plots present angular dependence 
of the singlet yield for different environmental noise levels 

[r = o(«),r = io 4 («),r = io 5 (A),r = 4 x io 5 (*),r = 

8xl0 5 (T), andT = 10 7 (*)]. On the right plots we present the 
magnetic sensitivity defined in Eq. ([7]). They reveal resonant 
character of the sensitivity for noise levels close to T = 4 x 10 , 
which correspond to optimal decoherence time of the order of 
a microsecond. Moreover, the upper right plot shows that the 
sensitivity is increased in the presence of environment. 



sensitivity. 

In summary, we have used the results of two differ- 
ent behavioral tests preformed with European Robins to 
estimate the average life time of the radical pair tak- 
ing part in the avian magnetoreception. Unlike a recent 
study which took into consideration the result of only 
one behavioral test and estimated the average life time 
close to 100 [is [l5j, our estimation of the life time is 
about few microseconds which agrees well with experi- 
ments. As the most important result of the present work 
we consider identification of a parameter regime where 
the presence of environment enhances performance of the 
chemical compass. Similar enhancement caused by envi- 
ronment is also found in studies of energy transfer during 
photosynthesis I6|.l35lj , and very recently in avian magne- 
toreception [34], |36| . These suggest that Nature might be 
optimizing performance of some biological processes by 
utilizing inevitable noise present in environment. More 
insight into this conjecture can be obtained from further 
studies of identified here resonance that magnetic sensi- 
tivity displays as a function of environmental noise. 
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the noise level T. Our aim is to study influence of the 
noise level on the magnetic sensitivity of the compass, 
defined as 

A s = $£ ax - $f n , (7) 

where $™ ax (<£>g ln ) is the maximum (minimum) singlet 
yield optimized over angle 9. 

Fig. [5] presents results of calculations of the singlet 
yield and the magnetic sensitivity for different environ- 
mental noise levels. For noise levels below V 10 4 we 
find that the sensitivity is practically independent of en- 
vironment. Therefore, for decoherence time To = s» 
100 [is and longer, the sensitivity is practically the same 
as if there is no environment. For shorter decoherence 
times (bigger T) the noise does influence the sensitiv- 
ity and for the HF strength of a = 1.0 [is we find quite 
counterintuitively that the sensitivity in the presence of 
noise is better than without noise (r = 0). This is not 
so for a — 0.5 [is, but in both cases the sensitivity is 
displaying resonant behavior as a function of T with its 
local maximum at the noise level r = 4 x 10 5 , see sup- 
plementary material [3l| for further evidence that this 
resonant noise level does not depend on k. Therefore the 
best magnetic sensitivity for a = 1.0 [is case takes place 
when the decoherence time To = = 2.5 [is. Here 
longer decoherence time is not useful as it can spoil the 



Appendix A 

In the main text, we use two behavioral experiments to 
estimate the life-time and the coherence time of the rad- 
ical pair (RP). We set the static magnetic field at 47/iT, 
which is the local geomagnetic field at Frankfurt, the site 
of all behavioral experiments. The life-time of RP is esti- 
mated using resonant RF field with strength 150nT, the 
RF field which was used in a recent behavioral experi- 
ment. In the same experiment, birds' behavior was also 
studied for other RF field strengths and also for stronger 
static magnetic field. This supplementary material (SM) 
discusses results which we obtain using some of those dif- 
ferent parameter values. They show that our estimations 
of the time scales are robust and do not depend much 
on the strength of RF field, and our identification of the 
resonance of magnetic sensitivity does not depend much 
on the life time of the radical pair. 

1. Life-time estimation for RF field of 470nT 

In Fig. [3j we show the angular dependence of the sin- 
glet yield for different decay parameter (k) values [l3| 
and for three different kinds of magnetic field environ- 
ment: in the geomagnetic field of 47/xT, in 30% weaker 
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FIG. 3. (Color online) Life-time estimation of a radical pair 
using 470nT RF field. Each plot presents angular depen- 
dency of the singlet yield in presence of local geomagnetic 
field of 47 fiT (•), in 30% weaker field (■), and for the geomag- 
netic field augmented with the resonant 470nT RF field (a). 
Calculations are performed for two strengths of the hyper- 
fine coupling: slightly bigger than the geomagnetic strength 
(a = 1.0 /is, three lower plots) and slight smaller than the ge- 
omagnetic strength (a = 1.0 [is, three upper plots). On plots 
from left to right we display the growing value of the recom- 
bination rate k used in numerics. The radical pair life-time is 
the inverse of k, and can be estimated using the argument of 
the main text. 



FIG. 4. (Color online) Magnetic sensitivity of avian compass 
in presence of environmental noise. The upper plots present 
results of calculations for hyperfine coupling strength of a = 
1 us and k = 2 x 10 5 , the lower plots are for a — 0.5 /is and 
k — 1.5 x 10 5 . The left plots present angular dependence 
of the singlet yield for different environmental noise levels 
[F = 0(«),F = 10 4 (B),F = 10 5 (A),r = 4 x io 5 (*),r = 
8x 10 5 (T), and F = 10 7 (*)]. On the right plots we present the 
magnetic sensitivity defined in Eq. (7) of the main text. They 
reveal resonant character of the sensitivity for noise levels 
close to r = 4 x 10 5 , which correspond to optimal decoherence 
time of the order of a microsecond. Moreover, the upper right 
plot shows that the sensitivity is increased in the presence of 
environment. 



magnetic field, and in the geomagnetic field augmented 
by weak resonant RF field of £? r f = 470nT (1% strength 
of the local geomagnetic field). Applying the argument 
of the main text gives estimation of the life-time 1.67/is 
for a — lfis, and 2[is for a = 0.5^s. 

2. Resonance of magnetic sensitivity for life-time of 
about microsecond 

Fig. 2] shows qualitatively similar influence of noise on 
the sensitivity of the avian compass, as that of Fig. 2 
of the main text (values of k are different here). More- 
over, the maxima of the resonance for all of these plots 
correspond to the same noise level. 

3. Life-time estimation for RF field of 47 nT 

We have also estimated the life-time of RP using much 
weaker resonant RF field of strength f? r f = A7nT (0.1% 
strength of the local geomagnetic field). Fig. [5] shows 
that for this case, estimated life-time is ll/us for a = l/j,s 
and 12.5/Lts for a = 0.5/is. 



4. Life-time estimation for stronger static field 

In Ref . [lH , birds were also tested by exposing them to 
artificially created magnetic field with strength twice as 
big as the local geomagnetic field. Birds adapted them- 
selves in this stronger magnetic field environment, and 
hence the artificially created stronger magnetic field is 
now the effective local geomagnetic field for the birds. 
The experiment shows that, in this situation, birds are 
again disoriented by applying RF field with frequency 
resonant with the Larmor frequency corresponding to 
this stronger magnetic field. In addition, according to 
Ref. [2l|, the 30% functional window should also exist 
for this case. Following these, we have performed similar 
analysis for the stronger 94^T artificial magnetic field to 
estimate the life-time. For this purpose, we use resonant 
RF field with B ri = 150nT, the same RF field which is 
used in main text. Fig. [5] presents results of these calcu- 
lations which reveal the life-time of 2.5/xs for a — 1/xs and 
5/xs for a = 0.5/is. Both values are again close to those 
estimated in the main text and show their robustness. 
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FIG. 5. (Color online) Life-time estimation of a radical pair 
using 47nT RF field. Each plot presents angular dependency 
of the singlet yield in presence of local geomagnetic field of 
47 fiT (•), in 30% weaker field (■), and for the geomag- 
netic field augmented with the resonant 47 nT RF field (a). 
Here again calculations are performed for two strengths of 
the hyperfine coupling: slightly smaller than the geomagnetic 
strength (a = 0.5 fj,s, three lower plots) and slight bigger than 
the geomagnetic strength (a = 1.0 ps, three upper plots). On 
plots from left to right we display the growing value of the 
recombination rate k used in numerics. Here we estimate 
the average life-time 11/xs for a = 1.0 /is case and 12.5/xs for 
a — 0.5 fis case. 



FIG. 6. (Color online) Life-time estimation of a radical pair 
for the static field B = 94/iT using 150nT RF field. Each 
plot presents angular dependency of the singlet yield in pres- 
ence of artificial magnetic field of 94fiT (•), in 30% weaker 
field Bo = 65.8/iT(B), and for the artificial magnetic field 
augmented with the resonant 150 nT RF field (a). As usual, 
calculations are performed for two strengths of the hyper- 
fine coupling: slightly smaller than the geomagnetic strength 
(a = 0.5 fis, three lower plots) and slight bigger than the ge- 
omagnetic strength (a = 1.0 /is, three upper plots). On plots 
from left to right we display the growing value of the recom- 
bination rate k used in numerics. We estimate the average 
life-time 2.5/xs for a = 1.0 /j,s case and 5fis for a = 0.5 /xs 
case. 
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